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LINEARIZED THEORY OF STAGNATION POINT HEAT AND MASS 


TRANSFER AT HYPERSONIC SPEEDS 
By Kenneth K. Yoshikawa 
Ames Research Center 


SUMMARY 


An explicit closed-form solution has been obtained for the compressible, 
•viscous shock-layer energy equation with mass addition or suction. The solu- 
tion is applicable at very high speeds, when the interactions between the 
heated gas behind the shock wave and the layer of injected or ablated gases 
from the surface must be considered as well as the effect of heat conduction 
across the shock layer. The solution consists of the energy equation 
decoupled from the other conservation equations and then linearized by the 
introduction of constant -values for the thermal and transport properties . 

Results of using the solution show explicitly that the effect of gas 
injection on convective heat transfer depends on the Stanton number without 
mass addition and on the relative molecular weights, and thermal and transport 
properties of the free-stream and injected gases. The solution is 
substantiated by comparisons with more exact solutions. 


INTRODUCTION 

A number of investigators (e.g. , ref. 1) have suggested that reducing 
the bluntness of a vehicle will significantly decrease the severe radiative 
heating it encounters as it enters the earth's atmosphere at very high 
speeds. However, reducing the effective bluntness increases the convective 
heating to such a body. Thus, there is renewed interest in the problem of 
convective heat transfer, including the effects of mass addition employed for 
cooling. 

Theoretical investigations of the effects of mass addition on convective 
heat transfer have been presented in references 2 through 20. The injection 
of nonreacting gases is generally treated by one of two fundamentally differ- 
ent approaches; (l) boundary -layer theory (refs. 2-5) and (2) shock-layer 
theory (refs. 6-8). 

At altitudes and velocities associated -with high-speed entry, it is 
necessary to employ shock-layer theory to analyze properly the problem of 
convective heating with mass addition. For high Reynolds number flow, Howe 
and Sheaffer (ref. 6) pointed out that, contrary to the results obtained from 
boundary -layer theory, the heat reduction due to mass addition was not simply 
correlated as a function of mass addition, rate. For low Reynolds number 


flow^ Chen et al. (ref. 8) pointed out that the heat reduction due to mass 
addition was considerably different from that obtained with boundary -layer 
theory. Although the investigations based on shock-layer theory are useful^ 
each problem^ because of its complexity, requires a separate numerical solution. 

With regard to the heat transfer, much of the complexity associated with 
the solutions to the shock-layer equations can be avoided if the energy equa- 
tion is decoupled from the continuity and momentum equations . This approach 
seems plausible in light of recent results (refs. 21-23). Hanley and Korkan 
(refs. 21 and 22) showed that for inviscid flow, the normalized mass-flow 
distribution along the stagnation streamline was insensitive to radiative heat 
transfer in the shock layer. Matting (ref. 23) showed that for viscous flow 
along a cold wall, the mass-flow gradient in the boundary layer was essentially 
constant at the corresponding stagnation-point value. 

The purpose of this report is to present a method for solving the shock- 
layer energy equation, decoupled from the continuity and momentum equations, 
by means of the results from references 21 to 23. (A similar method involving 
only radiative heating was presented in ref. 2k.) Closed-form solutions of 
the energy equation are presented to demonstrate the dependence of convective 
heating on the appropriate physical parameters. In order to substantiate the 
present method, these solutions are compared, wherever possible, with numeri- 
cal solutions obtained from the coupled equations. The primary effect of heat 
conduction in the shock layer on convective heating at the stagnation point is 
also incorporated in the analysis. 


SYMBOLS 


a 

B 


D 

erf 

H 

h 

h 

j 


parameter defined by equation (lie) 

Xw 

blowing parameter, — 

characteristic blowing parameter (eqs. (B18) and (B19)) 
mass fraction 

specific heat with constant pressure 
diffusion coefficient 
error function 

1 2 

free -stream total enthalpy, h^ + — 
enthalpy 


dimensionless enthalpy, j 

total enthalpy, h + ^ (u^ + v^) 


2 



k 

kb 

kw 

L 

M 

Uu 

ijf 

Pr 

P 

q. 

R 

Re 

r 

St 

(St)^ 

T 

T 

u 

V 

V 
X 

y 

P 

r 


total theimal conductivity 
constant defined in equation (b 4) 
constant defined in equation (BIO) 

(x’)** 

effective ratio of mass-flow gradient (eq. (B9) ) ^ — 

(X^)* 

shock standoff distance 

molecular weight of cold gas (at wall temperature) 

Nusselt’s numher (eqs. (B25) and (B26) ) 

factor of proportionality defined in equations (l 91 >) and ( 29 b) 

Prandtl numher 

pressure 

heat flux 

body radius 

Reynolds number 

radial distance defined in figure 1 
Stanton niunber, 

Ir 

Stanton number for radiative heat flux, 

temperature 

T 
■^s 

X component of velocity 

X 

value of u at the shock, ^ 
velocity 

y component of velocity 

distance along body 

distance normal to wall 
p V L 

Peclet number, , (eq. ( 6 b)) 

modified factor of velocity gradient (eqs. (Bk) ) 
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boundary-layer thickness 


€ 

•n 

A 




p 



X 

K 


* 

** 


A,B 

c 

cr 

f 

m 

o 

r 

ref 

k 


characteristic thickness 

P 

density ratio across the shock wave, ^ 
effective heat of ablation 
normalized distance from wall, ^ 
parameter defined in equation (l 9 b) 
viscosity 
density 

scaling factor for heat transfer normalized by value for 

pv 


= 1 


Tg = 15 , 000 ° K, R = 1 ft 
dimensionless mass-flow parameter. 


Sx 


P V 
^00 00 


absolute value of 5 — near the wall 

orj 

convective heat reduction with gas injection. 


1 , 


■CO 


Superscripts 

Integrated quantity for the condition Tg » (e.g., eq. (9a)) 

quantity for the condition Tg T^ (eq. (B5)) 


Subscripts 

species A,B 

convective 

critical 

foreign gases or equivalent foreign gas 

mixture 

no blowing 

radiative or reference conditions 
reference properties at T = 15,000° K 


atm. 




s condition just tehind shock wave 

w wall 

6 houndary-layer edge 

<» free-stream gas or conditions 


ANALYSIS 


Conductive heat transfer will he considered with a simplified energy 
equation for the viscous , compressible flow in the neighhorhooa of a stagna- 
tion point. The thermodynamic and transport properties of the shock-layer 
gas will he introduced as averaged quantities. Solutions for the resulting 
linearized equation will he obtained for a body from which transpiration or 
ablation gases are emanating in hypervelocity flight. The calculation of the 
convective heat -transfer rate at the wall will he emphasized. 


Assumptions and Geometry 

The following assumptions are basic to the analysis (fig. 1 shows the 
geometry used in the analysis) . 

1. Steady axisymmetric stagnation -point flow 

2. ThermDd 3 mamic equilibrium 

3. No chemical reactions between free-stream and injection or 
ablation gases 

i^-. Thermal radiation uncoupled (convective heat can be derived 
separately by means of the driving enthalpy) 

5. Laminar continuimi flow 

6. Thin shock layer, L/R « 1 

7. Small blowing rates, p^v^ « 

8. Thermal diffusion and diffusion-thermal neglected 



Figure 1.- Flow geometry. 


5 



Basie Equations 


Conservation equations for steady flow of multicomponent chemically 
reacting gases in the shock layer have been derived previously (see, e.g., 
refs. 5^ 25 - 27 ). The resulting equations (discussed in detail in ref. 26) for 
the conditions of the present analysis are: 


Continuity equation 



^ (pur) + ^ (pvr) = 0 

Momentum equation 

6u , 6u ( ^ 

P^ 6x P^ Sy 

3 \ 6 / 6u\ 

c/ ^ 6y 5yj 


3 ^“ 0 

oy 


Energy equation 

B .1 . Bi B b7 

PU 5 ^ + PT ^ 

) ^ (‘Ir) 

Diffusion equation for the binsiry system 

6c|. 6cp ^ / 

Sr + 6y = Sy \ 

' be A 

By j 


(la) 


(lb) 


(ic) 


( 2 ) 


(3) 


Since the major portion of this study is concerned with convective heat 
transfer to the stagnation point, the energy equation for the shock-layer flow 
is further simplified, for pu(6j/6x) « pv (6j/6y) and (l/2)(u^ + v^) « h. 


6h _ a 
6y 6y 


/jl 

\cp bjj 



(4) 


Equation (4) can be decoupled from the continuity and momentum equations (l) 
by the introduction of a normalized mass-flow distribution which is 
Independent of the heat transfer. 


Normalized Mass-Flow Distribution 

The mass-flow distribution along the stagnation streamline must be 
determined to decouple the energy equation. For inviscid, nonadiabatic flow, 
Hanley (refs. 21 and 22) showed that the normalized mass-flow distribution, 

X = pv/p V , is relatively insensitive to radiation transport. For 
boundary -layer flow. Matting (ref. 23 ) demonstrated that the mass-flow 
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gradient through the highly compressed "boundary layer was essentially constant 
These tuo results suggest that the mass-flow profile throughout the shock 
layer (-viscous and nonadiabatic) is also relatively independent of the heat 
transfer^ an ass-umption substantiated by a number of published papers (refs. 6 
7, 265 28 , and 29 ). (Note that the flow-field solutions in refs. 6, 7 , and 
26 include mass addition in the stagnation region.) Thus, the effective mass- 
flow equation -will be written as 

pv 

X = Y - x^ - “ (1 “ X-^ + x-y)T) ( 5 ) 

l-'oo 00 

where ti = y/L is the distance from the wall normalized by the shock stand- 
off distance. It should be noted that X-Jf is the absolute value of (Sx/Stj) 
a short distance away from the wall since the wall -value is zero. One is 
justified in this approach because Xw does not directly affect the heat 
-transfer (see eq,. (8)) and 8 x/8t) is essentially constant through the bound- 
ary layer except very near the wall (see, e.g., refs. 23 and 26, figs. 8 
and 2, respectively) . It will be assumed that Xv independent of the 
blowing rate but a function of the body shape and density ratio across the 
normal shock. 


General Solution of Energy Equation 


From equations (4) and ( 5 ) the dimensionless energy eq,uation for the 
shock layer becomes 


dh 


1 dh 


n\ d(St). 


^ df) dri yP dt) J 

where dimensionless quantities are defined as 

h 


dT) 


h = 


H - h, 


w 


(St)^, = 


P = 


Pco^ooCH - K) 
PooVooL 


k/c. 


•X 


J 


(6a) 


(6b) 


and H is the adiabatic enthalpy behind the shock wave. For simplicity, 
equation (6a) may be solved by treating separately the effects of conduction 
and radiation on the flow enthalpy. The case where radiation is predominant 
has been studied in reference 24. The analysis in this report is concerned 
with the case where conduction is predominant. 


7 



The energy equation and boundary conditions, neglecting radiation become 


X 


^ ^ A. (l. ^ ) 

d-q dri \ ^ dT] / 


(7a) 


h = at T] = 0 


h = hg at T] = 1 


(7b) 


A general solution to equation (ja) and boundary conditions (tI)) may be 
written /* t] ^ 

_ _ r px dT] 

h - h^ = (St) / pe ° ^ dri^ 


(8a) 


where 


(St) = i (fl) . 


h - h 

s w 


pe 


n X 

Jo 


(8b) 


df] 


Equations (8) are nonlinear solutions since p depends upon the unlmown 
function h. Generally, the analytical evaluation of this equation is 
extremely difficult. 


Linearized Solution of Energy Equation 

A linearized solution to the energy equation can be given explicitly when 
the Peclet number p is replaced by an average value. For the condition of 
constant pressure and T^ « Tg, let 


and 



(9a) 


P*(Tg) 


P 


V L 
00 00 


(k/cp)* 


(9b) 
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The linearized solution thus becomes 



where 



(lOa) 


(10b) 


A vigorous mathematical justification of the substitution of this average 
function in the nonlinear energy equation is beyond the scope of this report. 
However, the method is substantiated later where the numerical results of con- 
vective heat transfer from the solutions of other investigators (refs. 6, 2.6, 
30 ) and also from the exact solution of equations (8) show excellent agree- 
ment with the present linearized solutions. Incorporating equation (5) in 
equations (lO), one obtains 


1^1 „ “P*[-Xvn+(i/2)Xw’i^+(i/3)(i-x4+Xw)'n^] 


(lla) 


when 3* >> 1, the entire exponential tends to zero. Furthermore, when t) 
is small, the last term in the exponent can be neglected. This allows the 
integral of equation (lla) to be expressed in terras of the error function 


_ (h-s 
h = 


hjerf (PX^)* T 1 - a + hg erf a + h^ erf (pX^)* “ a] 


erf 


[J2 




+ erf a 


(llh) 


where 

Furthermore, for X-^ small and p* » 1, 



(11c) 
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It follows then from equations (ll) that 


1 1 


“ Ji 1 + erf a 


(l2a) 


The convective heat is^ by definition^ 


^ (-L = Poo^oo^^ ~ 

\cp Sy/^ \dTi/^ 


(12b) 


Thus, 


2 


■(i/2)(p/X^)*X^ 


= Poo^J^S - l^w) 


1 + erf 


(4) 


’•w 


(13) 


For the no Injection case X-^ ^ 0, and equation (13) can be written 


4co = /I (r) pjJ\ - K'l 


(U) 


where 


P V L 

* _ '^00 00 o 


A blowing parameter and a Stanton number for the no injection case can be 
defined as 


B = 


X-y 


^CO 


^ P V (H - iv) 


f 3 T \ Pq / 


for :— « 1 
ks 


(15a) 

(15b) 


Equations for Heat Reduction 

The heat reduction due to foreign gas injection is obtained in terms of 
the blowing parameter B by dividing equation (l3) by equation (l4) and 
using the definition (l5a). 
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0 

11 

H 

1 L (’Ck/cp)* 

^0 (x>/cp)* 
e 

II 

JLo (XJk/cp)J 

P /l L 

^ Jk L, ()Ck/cp)- ® 


(16) 


where subscripts t» and m refer to properties of the free -stream gas in the 
shock layer for no injection and properties of the mixture of gases in the 
shock layer for injection^ respectively. The closed form solution (eq. (l6)) 
can be used to deduce the ratio of the parameter (X-^k/cp)*/(X-^/cp)* from 

either experimental data or the results of numerical calculations (e.g., 
refs . 3 • 

Identical-gas injection . - With the assumption that is independent of 

blowing, it follows directly for identical-gas injection that 


(^/^p)co 

«k/cp)^ 


Then i); is written 

- i ^ 

e ^ 

1 + erf ^ B 



(17) 


Equation (I7) provides remarkable agreement with more exact solutions as 
shown in the Results and Discussion section. 

For the case L/Lq 1 , such as for a large body with a thin boundary 
layer, equation (17) may be expanded in a series as 

t 1 - - B + . . . (18) 

It 

Equation (18) compares very well with the empirical formula suggested by 
others (refs. 2, 3, 9 > 10 ). 

Single foreign gas injection .- The property, (k/cp)^ in equation (I6) 
refers to the mixture of gases in the shock layer and as such it depends on 
both the individual thermal and transport properties and the local mass con- 
centrations of the species. Therefore, an exact evaluation of this property 
would entail a solution to the species equation and moreover necessitate the 
computation of high temperat'ure transport properties which requires physical 
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constants still not well established (refs. 11, l4, l8, 31"3T)* Since the 
primary interest here is to obtain a simple heat-transfer formula for engi- 
neering use, the above difficulties will be avoided by introducing the 
following approximation 


(X>/cp)^ 

(^w^/ <^p)j^ 


(xWcp)* 

(X^/cp)f 
_(x.-^/cp) j 


for 0 < B £ 


for B > Bj, 


where B* is the blowing parameter for the boundary- layer blow-off defined 
in appendix B. Eq^uation (l9a) has the proper physical behavior since it 
weights the properties of the foreign gas according to the magnitude of B 
for B ^ B* and is constant for B > B* where the foreign gas is dominant 
over the main stream gas. 

Since (k/c^)* , the parameter A may be defined by 


(x^/c ) 




(X^/Cp)f s/^ 


(19b) 


where Nf is a weighting factor which depends on the type of foreign gas . 
By the present method in conjunction with the numerical results of refer- 
ence 3^ bhe factor Nf was found to be (for air as the main stream) 




for diatomic and polyatomic gases 


(19c) 


Derivations of these values will be discussed in detail in the Results and 
Discussion section. 

Equation (l6), in conjunction with equations (l9a) and (l9b), therefore, 
leads to 
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Ilf 


fc 


(B/B*)' 


A 


-1/2 


_ i JL 7 ,(b/b*) ^2 






t 


fe^)' 




— — 

1/2 

J 


1 L 



1 + erf 

7 T ^ 

Jt Lq 

B 


1 L 

•N T^2 



L 

XB 


1/2 

0 

e 




for B < 


(20a) 




for B > B* 


(2OB) 


Equations (20) provides excellent agreement with the results of reference 4. 


It is interesting that if the empirical form of the blowing parameter 
suggested by a n-umber of authors for the case of foreign gas injection (i.e., 
replace B by (l^/Mj.) B) is used, equation (I7) immediately leads to 



where is the molecular weight of the foreign gas. Equation (2I) 

provides very satisfactory results for most single foreign gas injectants 
(hydrogen excepted) . 


Mixed foreign gas injection. - When injected gases consist of two or more 
gases, the calculation of the transport properties is further complicated. 

A reciprocal formula for the equivalent transport properties is applied, 
that is , 


1 ^B 

(X-^k/cp)^ (Xvk/cp)j^ (X-j^k/cp)-g 


(22a) 


where A and B refer to the Injected gases and are their respective mass 
fractions. Incorporating equation (22a) in equation (l9b). 


X = 



^B 


Mb %/ 


(22b) 
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This may he ■written in terms of the equi'valent molecular ■weight and N^. 
as 

A=^Nf (22c) 

where 


3^ ^ fA ^ 

% Ma % 



Substituting equations (22) into equations (20)^ >lr for the case of the mixed 
gas injection is obtained. As can be recognized from equations (22 ) , the 
equi'valent molecular weight consisting of two or more different gases does 
not itself pro^vlde the same result for a single gas with the same molecular 
weight. It will be shown later that equations (20) and (22) pro'vade close 
agreement with the mmierlcal solutions of reference b when more than one gas 
is injected. 

Ablation . - For forced injection the blowing rate pareimeter is given a 
priori and the solution is straightforward. For ablation, however, the blow- 
ing rate parameter depends on the heat transfer to the wall, and the solution 
is, therefore, coupled. The mass loss rate of ablation vapors per unit area 
is given by 


dm 

dt Pw'^w 



(23a) 


where ^ is the effective heat of ablation. With the definition of B 
given in equations (l5), equation (23a) may be rewritten 






- K) 


B 


( 2313 ) 


Ideally, equations (l6) and (23I3) may be combined to eliminate B (or i]/) and 
to obtain ijr (or B) for ablation. However, a more practical approach is to 
obtain graphical solutions by superimposing equation (23!)), which represents 
a family of straight lines with the slope ?/Ch " hy) , on equations (20) for 
ijr calculated for the forced injection. 


Effect of Heat Conduction Behind the Shock 

The enthalpy immediately behind a strong shock wave will be significantly 
altered by heat transfer at the ■wall when conductive heat becomes a major 
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factor in the energy balance ac 2 X)ss the bow shock wave. The effect of 
conductive heat on the energy balance is to reduce the flow enthalpy behind 
the shock; this phenomenon is designated here as the "postcoollng effect," 
in contrast to the preheating effect which Increases the flow enthalpy behind 
the shock (ref. 24). Only preliminary effects of heat transfer on the vis- 
cous flow field and on the convective heat at the stagnation point will be 
discussed (see Results and Discussion). A brief ainalysis of the reduction in 
heat transfer due to the postcooling effect is described in appendix A. 


Approximate Equations for the Flow Properties 

The various physical quantities required to evaluate the linearized 
solutions are given, or derived, in detail in appendix B: These quantities 

are the mass-flow gradient, shock-layer thickness, boundary -layer thickness, 
boundary -layer blow-off parameter, and the heat -transfer parameter. 


RESULTS AM) DISCUSSION 


Calculations of the average transport properties, the convective heat- 
transfer rate at the stagnation point with suction or foreign gas injection, 
the heat -transfer parameter in the vicinity of the stagnation point, and the 
mass-flow gradient for the higher wall temperature are presented in this 
section for atmospheric flight conditions at very high speeds. Also con- 
sidered is the postcooling effect on the convective heating rate at low 
Reynolds numbers. 


Thermal and Transport Properties 

The thermal and transport properties for air reported by Hansen (ref. 3I) 
are employed because the present analysis will be compared with numerical 
solutions using these properties, and changes in the properties at the higher 
temperatures as reported recently (refs. 34-3?) should not significantly 
affect the conclusions of the present analysis. 

Normalized values of k/cp are shown in figure 2(a) as functions of 
temperature for various pressure levels. (For simplicity in plotting, the 
reference condition is based on T^^p = 15,000° K.) Also shown are the 
averaged values, (k/cp)*, calculated by equation (9a-) • The parameter (k/cp)* 
can be approximated by a straight line which passes through the origin and a 
reference value at a lower temperature. (See tables I and II for mmerical 
value s . ) 

Figure 2(b) presents the average parameter (k/cp)* plotted against 
enthalpy at 1000° K intervals. It can be seen that (k/cp)* may be approxi- 
mated by the square root of the enthalpy. This also follows from the previ- 
ous figure since temperature is approximately proportional to the square root 
of the enthalpy (see, e.g., fig. 2 of ref. 24). 
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Figure 2 . - Thermal and transport property 
parameter for air. 



(h) Average -property variation with enthalpy. 
Figure 2.- Concluded. 


The results found from figure 2 
hear an important consequence; namely, 
the thermal-transport property param- 
eter can he approximated hy 



(24a) 


where hj. is the reference enthalpy 
at a lower temperature (a constant for 
all gases). ^ Incorporating equa- 
tion (24a) with equation (l4) , one can 
derive the important equation (see 
eq. (B22)). 







^w) 


(24h) 


Equation (24h) shows that the 
convective heat transfer at high 


^That is, ^ 1700 Btu/lh (=“ 4 Mj/kg) . However, the reference enthalpy 
at a higher temperature which provides the same result in equations (24) is 
hj;. 27^000 Btu/lh (=“ 63 MJ/kg) for all gases. As will he recognized later, 
the parameter used in equation (l9h) can he deduced more precisely from the 
figures at the higher h^ than the lower h^. 
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temperatures can be expressed in terms of gas properties at a lower tempera- 
ture. Marvin and Pope arrived at the same result for air and other gases. 
Equations (24) further provide the following important relations: Equa- 

tion ( 19 b), with equations (24), becomes 


[ (X>/ Cp ) * ] ^ ( qco Je/Ps) ^ 

■ ■ ■ ■ I I 

[ (^^jl^/cp) f ( Ico •J E/Pg ) f 


(25a) 


Since (k/c.^)* CC Jii at the lower temperature, where M _is the molecular 
weight of the cold gas, the factor of proportionality Nf is related by 


^ (lco>s/E/Ps)| 


( 25 b) 


The values of Nf given in equation (l9c), and equation (22c) for gas 
mixtures, can therefore be determined from figure 5 in reference 3 . 


Results Without Postcooling 


Laminar convective heat transfer without in,1ection .- The results for no 
postcooling obtained from equations (B3), (.B 6 ), and CB21) are shown in fig- 
ure 3 (a) along with the predictions of Hoshizaki (ref. 30 ), Fay and Riddell 
(ref. 38 ), and Howe and Viegas (ref. 26). Typical results from other studies 
may be found in references 39 to 43 . The present results show a slight but 
noticeable dependency on pressure level. They agree well with the numerical 
results of Howe and Viegas for Pg = 1 atm, and approach the result of 
Hoshizaki for lower pressure. 


The heat -transfer parameter Nu/^Rej calculated from equation (B2T) 
assuming (X-^)* — I /2 is shown in figure 3(b) for Pg = 1 atm. The present 



(a) Convecti-ve heating rate. 


Nu 

//Re X 



O PRESENT METHOD, = lOtm 

T^= 1.500" K 


8 12 
Va>, km/sec 


(b) Nusselt nmnber. 


Figure 3.- Stagnation-point heat transfer for air. 
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analysis provides remarkably good agreement with the results of Fay and 
Riddell for the lower flight speeds and with the results of Howe and Viegas 
for the higher flight speeds. 


Mass -flow gradients . - The effect 
mass-flow gradient is demonstrated in 
The quantity k^^ the effective ratio 



Figure it-. - The effect of the wall ten^ierature on 
mass flow gradient (air) . 


of wall temperature on the normalized 
figure 4 for specified conditions, 
of the mass -flow gradient for the hot 
wall to that for the cold wall^ 
deduced from reference 23^ is compared 
with the present result (eq. (b 8)) for 
the condition of the lower stagnation 
temperature. The effective mass-flow 
gradient decreases as the wall tem- 
perature increases because of the 
effect of higher viscosity in the 
boundary layer. Also presented is 
the result for the higher stagnation 
temperature (Tg = 15 ^ 000 ° K, 

Pg = 1 atm) which demonstrates a 
further decrease in the effective 
velocity gradient. 


Heat reduction - identical gas in.jectlons . - A comparison of the present 
solution (eq. ( 17 )) for p^H" = 1 with the more exact numerical solutions of 
Howe and Sheaffer (ref. 6), and with the correlation of numerical solutions 



Figure 5 . - Heat reduction with identical gas 
injection. 


by Marvin and Pope (ref. 3 ) ^ whose 
results are identical to those of ref- 
erence 2, is presented in figure 5 
for identical gas injections. The 
present results show excellent agree- 
ment with all of the numerical solu- 
tions. Also presented is the linear 
equation (18) which agrees with the 
empirical equations of references 9 
and 10 and is very close to the 
linearized solution for B < 1. Note 
that the ilr function is independent 
of the boundary temperatures Tg and 
(ei- ( 17 ))- However, the present 
solution shows that the ilr function 
is not a single curve, but depends 
upon the stagnation-point Stanton num- 
ber (st) or the product of the 

ambient density and body radius p R 
as shown next. 


Scaling factor for correlation . - As the product of the density and the 
body radius decreases, the mass transfer becomes more effect ive in reducing 
the heat transfer, as illustrated in figure 6. The product p R in the 
figure has been normalized by the value for Pg = 1 atm, Tg = ^5,000° K, and 
R — 30 cm (1 ft) . The present solutions and those of Howe and Sheaffer 


18 




(ref. 6)^ agree well. A single 
correlation curve for the shock-layer 
theory c an he obtained only if the 
product p^R remains constant. This 
scaling effect on heat transfer is 
easily recognized from equations ( 17 ) , 
(B11) , and (B12) . Boundary -layer 
theory, on the other hand, provides 
one correlation curve of the blowing 
effectiveness, independent of the 
product p^R (e.g., refs. 2-4). For 
large values of p^R, the shock-layer 
and boundary-layer theories seem to 
provide similar results . The present 

results and those of Howe and Sheaffer and of Marvin and Pope (also refs. 2 and 4) 
agree well, as illustrated in the previous figure. This is a condition for 
which the shock-layer thickness changes only a small amount as a result of gas 
injection. Changes in the shock-layer thickness res ult ing from transpiration 
or gas injection are significant for low values of p^R (but still in the 
continuum regime without the postcooling effect) and are accompanied by 
changes in the velocity gradient at the stagnation point. The present shock- 
layer solutions account for the changes in the shock layer thickness, as do 
those of Howe and Sheaffer and of Goldberg and Scala. 


Figure 6.- Heat reduction for air with scaling 
factor. 



Figure 7.- Heat reduction in air with single 
foreign gas injection. 


Heat reduction - foreign gas 
injection . - Injection of gases dif- 
ferent from the free -stream gas 
affects the heat transfer in the man- 
ner shown in fig ure 7* The present 
results for p^R = 1.0 using equa- 
tions (20) are compared with the 
boundary -layer solutions from ref- 
erences 4 and 20. Very good agree- 
ment is obtained with the solutions of 
Anfimov and Al’tov for "single foreign 
gases" (non-mixed) over the entire 
blowing rate range (including suction) 
and equally good agreement is obtained 
with the solutions of Hoshlzaki and 
Lasher (ref. 20) for nylon phenolic 


ablation. The Injection of light gases is more effective in reducing the 
heat transfer provided no combustion takes place. 


The dependence of the effectiveness of transpiration upon the Stanton 
number (or free -stream density and body radius) is also a function of the 
molecular weight of the injectant gas relative to that of the free-stream gas. 


2Those numerical solutions that were misplotted in figure 5 of ref- 
erence 6 have been corrected in figure 6 of this report. 
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Figure 8.- Heat reduction for different free- 
stream gases. 



Figure 9." Heat reduction with mixed foreign gas 
injection. 


For example, figure 8 indicates that 
changes in the Stanton number affect 
the heat transfer more for the injec- 
tion of helium into air than for the 
injection of air into helium. Again, 
this is to he expected on the basis of 
changes in the shock-layer thickness; 
a gi-ven amotint of helium injected into 
an air stream increases the shock- 
layer thickness no re than the same 
aimDunt of air injected into a helium 
stream (see eq. (biO)). 

Figure ^ shows how injecting a 
"gas -mixture " affects heat transfer. 
For transpiration of a mixture of a 
light gas and a heavy gas, the gas 
near the wall, which is the heavy gas, 
predominates so that the mixture is 
less effective in reducing the heat 
transfer than a single gas^ of the 
same molecular weight. The dashed 
curve represents the result obtained 
for a single polyatomic gas of the 
same molecular weight. The results of 
the present method, using the recipro- 
cal of the eq^uivalent transport rela- 
tion (see eq.s. (22)), agree well with 
the numerical results of Anfimov and 
Al’tov. Anfimov and Al’tov proposed 
the simple formula: 


^ ^ J^l'^i (26) 

that is, the net heat reduction by the mixture is the cross product of the 
individual heat reductions for each species (ref. 4). However, this eq_uation 
is good only for small values of the blowing parameter. 


Effect of Postcooling on Convective Heat Transfer 

The effect of postcooling on the stagnation-point flow is to reduce the 
driving enthalpy, and thereby lower the convective heat at the body surface. 
(The vorticity at a low Reynolds number increases the convective heating 
(e.g., refs. 6 and 44), and thus tends to counteract the postcooling effect.) 
An exploratory study of this effect is presented here. 


^Single monatomic (polyatomic) gas if the mixture is made up with more 
than two monatomic (polyatomic) gases. 
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(Td) The enthalpy hehlnd the nonnal shock wave 
with air injection. 

Figure 10.- The effect of postcooling on heat 
transfer and the flow enthalpy (air: 

Pg = 1 atm, Tg = 15 , 000 ° K, and T^ =- 0° K) . 



Figure 11.- The effect of Stanton number on heat 
reduction, air. 


Stanton nimibers calculated by the 
exact solution (eq.s. (8)) and by the 
linearized solution of equations (lO) 
are presented in figure 10 (a) for the 
Initial condition of Pg = 1 atm, 

Tg = 15 , 000 ° K, =“ 0 , and the range 
of p* from 10“^ to 10 3 
(r =- 10"^ ~ 10^ cm) . The two calcula- 
tions are very near-ly identical. Note 
that the boundary condition behind the 
shock wave, hg, for smaller values of 
p*(hg) is now coupled with the post- 
cooling effect. The result of the 
simple formula without the post coo ling 
effect (eq. (l5b) with hg fn H) is 
also presented in the same figure, and 
agrees well with the exact solution 
for p* » 1. As can be seen from the 
figure, the effect of the post coo ling 
on the convective heating becomes 
important for the smaller body radius 
(or smaller Reynolds number) and at 
the limit for p* -» 0, Stanton number 
approaches unity, as for free molecule 
flow (or for a pointed nose) . 

Figure 10(b)_shows the variation 
of the enthalpy hg with post coo ling 
and injection obtained from the lin- 
earized solution (eqs. (lO)). The 
effect of the postcooling on the 
enthalpy becomes significant beyond 
the critical Stanton number, 

(St)cr ~ 0 . 32 . For a given value of 
(st)Q and increasing injection rates, 
the enthalpy is less affected by the 
postcooling since the conductive heat 
in the shock layer decreases consider- 
ably. After a sufficiently large 
blowing rate B ^ B* =“ 1.6, the 
enthalpy ratio is not significantly 
changed. 

Since the enthalpy varies with 
Stanton number and blowing rate, the 
heat reduction is altered in the man- 
ner shown in figure 11. The heat 
reduction, obtained from equations (lO), 
becomes less effective as the Stanton 
number increases beyond the critical 
number, (St),,!. O. 32 , and ilr 
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approaches the free molecule flow limit of iinity. Similar results were 
reported in reference 8. This same effect of the Stanton number on convective 
heat can he expected for foreign gas injections. 


CONCLUDING EEMAEKS 


Closed-form solutions have heen obtained which describe the behavior of 
the convective heat transfer to the stagnation point of a body with mass 
transfer. The .shock-layer energy equation was decoupled from the momentum 
and continuity equations by assuming that the distribution of mass flow (nor- 
malized by the free-stream mass flow) in the shock layer was independent of 
the heat transfer. The energy equation was linearized by the introduction of 
a constant thermal and transport properties for the free-stream and injection 
gases (i.e., it was found that convective heating can be simply calculated by 
an average Peclet number obtained by integration over the temperature range 
in the shock layer) . For the case of no radiation in the shock layer, the 
solutions indicate: 

1. The reduction of heat transfer by mass addition can be significant. 
The degree of heat reduction depends upon the mass addition rate and the 
Stanton number for no mass addition. For a given mass addition rate the heat 
reduction Increases with increasing Stanton number until the Stanton number 
reaches a value of about 0-5? at which point further increases in Stanton 
number decrease the heat reduction. 

2. The injection of gases of low molecular weight is most effective in 
reducing heat transfer. 

3- A mixture of injection gases is less effective in reducing heat 
transfer than a single gas with an equivalent molecular weight. 

The present linearized theory is well substantiated by comparisons with 
more exact solutions and provides explicit and simple correlation equations. 
The linearized theory of convective heat transfer, incorporating the lin- 
earized theory of radiative heat transfer (ref. 24), now provides the basic 
means for solving the coupling effect between convective and radiative trans- 
fer. The present method of analysis can be used to evaluate both the net mass 
transfer rate (including the diffusion which will be Important for lighter 
gases) and the shearing stress at the wall. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035^ Apr. 25, 1969 
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APPENDIX A 


THE EFFECT OF BOUNDARY CONDITIONS ON HEAT mAPTSFER (POSTCOOLING) 


The enthalpy immediately Behind a strong shock wave will he significantly 
altered hy heat transfer at the wall when conductive heat Becomes a major 
factor in the energy Balance across the Bow shock wave. Only preliminary 
effects of heat transfer on the viscous flow field and on the convective heat 
at the stagnation point will Be described in this appendix. 


Strong Normal Shock Relation Including Heat Conduction 

For simplicity it is assumed that the strong normal shock relation 
applies for the case of a conducting fluid. From a conservation of energy on 
each side of the shock wave 


pV (h +^V^')+q=pV (hs+^v2)+CL^ (A1) 

^00 00 \ 00 2 ^ 00 \ o 2 S / ' 


Since the conductive heat ahead of the shock is small (a =* O) and Vq « V , 
it follows that the enthalpy Behind the shock wave (including the Blowing 
effect) Becomes 


hs = H - Gs/p^V^ (A2) 

where H represents the adiahatlc enthalpy without conduction 

H = h + J (A3) 

Dividing equation (A2) By (H - h-^^) and differentiating equation (8) with 
respect to t), one obtains the Boundary enthalpy 


hg = 1 - (St)e' 


Jo 


(h„ 0) 


Also, from equation (8), hg Becomes 



and one can obtain the Stanton n\amBer from equations (a 4) and (A5) , 


(A4) 


(A5) 
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(St) = 


1 - Iv 



Jo 

pe d.T] + e 


(a6) 


The enthalpy behind the shock (eq^. (A5)) depends upon the blowing rate 
as well as the Stanton number. Equation (a 6) provides a closed form solution 
for p — Ps, a constant. When p = p*(hg)^ equation (a 6) can be solved by 
straightforward iteration. However^ th e d ifference in Stanton number due to 
changing p* is minor since p* — Pg* X/hJ for hg i‘ 0. 


Equation for Heat Reduction - Postcooling Effect 


The reduction in convective heat with gas injection will be signifcantly 
affected by the postcooling phenomenon. From equation (a 6) the heat reduc- 
tion is calculated by 


ilr = 


(St) 

(St)o 


(A7) 


It is observed from equations (A1) through (A7) that the major change in 
convective heat transfer due to mass addition for larger Stanton number (in 
effect, smaller Reynolds numbers) depends strongly upon the reduced enthalpy 
(driving enthalpy) accompanied by this postcooling process behind the normal 
shock wave. 
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APPENDIX B 


APPROXIMATE EQUATIONS FOR THE FLOW PROPERTIES 

The various physical quantities required to evaluate the linearized 
solutions are as follows: 


Mass Flow Gradient for Tg » T-^^ =“ 0 

According to the potential flow and the Newtonian flow theories, the 
velocity gradient at the stagnation point of an axlsymmetric "body is 
approximated hy (ref. 38 ) 


(Bl) 


For compressible flow and for Tg » T-^ it is found (ref* 23) that 


8y 



(B2) 


It follows that 





(B3) 


where 


V. 




eVoo 

(Lq/H) 

e 


1.187 

1 + 0.225 n/T 
1 









(b4) 


for a spherical body (see ref. 45 j see also ref. 46 for nonspherlcal bodies); 
7 is a modified factor introduced in reference 4T for velocity gradient. 


Mass-Flow Gradient (x^)** for Tg Ty 

When the wall temperature is comparable to the boundary -layer- edge 
temperature, the effective mass-flow gradient at the stagnation point 
decreases slgnficantly because of the higher viscous effect in the boundary 
layer . 


25 



Let the effective mass gradient at the wall for the case where Ts 
he (Xy)**^ the thermal and transport property parameter he 


^w 



(B5) 


Then^ the convective heat at the stagnation point for the no injection case 
is exactly analogous to that for Tg » Ty. It follows^ hy analogy with 
equation ( IL ) , that 


<lco = 



" ^w) 


(B6) 


It has been demonstrated (see refs. 26 , 30 , 38, and fig. 3(a)) that the heat 
transfer depends mainly on the enthalpy difference (hg - hy) ; therefore, it 
can he assumed, from equations (l^J-) and (b 6 ) , that 



(B7) 


Equation (B5) can he expressed 

(k/cp)* 

" (k/Cp)** 


in terms of (k/cp)*, and the ratio 

_ B - ^w 

- Ty[ (k/cp)*/(k/cp)*] 


(B8) 


where 



Tg 


From equation (B 7 ), the effective mass -flow gradient for general wall condi- 
tions becomes 


(x;)** “ i%(xy* (B9) 

This approximation is shown to he valid (fig. 4). 

Shock -Layer Thickness 

The shock-layer thickness L depends upon the mass injection rate. In 
adiabatic flow it is correlated hy a simple expression (refs. 46 and 48) 
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(BIO) 


where is the shape factor that depends upon the hody shape; for example^ 

reference 1+8 gives kb — 1.0 for an axisyrametric stagnation point. Using 
the definition for (St)Q given hy equation (l5a)^ equation (BIO) can he 
rewritten as 


+ (BID 

where the Stanton number at the stagnation point for the no injection case can 
he expressed simply as (e.g., see eqs. (ijb) and (B22)) 

(St)o cc (B12) 

Note that in this analysis the heat reduction equation for \|/ {see, e.g., 
eqs. ( 17 )) provides a correlation curve only if the Stanton number (St)Q or 
the product (p^oR) remains constant. 


Boundary -Layer Thi ckne s s 

An important characteristic dimension in the viscous layer is the ther- 
mal boundary -layer thickness 6. In this report the thermal thickness is 
arbitrarily defined as the distance from the wall where the driving enthalpy 
is 99 percent of its original value, that is, for y = S, 


^5 ~ ^ 
^s “ ^ 


0.99 


he 


(BI 3 ) 


Associating equation 


(BI 3 ) with equation (lib), 


one obtains the result 


6 

L 



^w 

/ • \* 

(k v) 


for « 1 


(B14) 


where the thermal layer thickness for the no injection case incorporating 
equation (l7b) is 



1.82 ^/^ 

I .* 

(kw) 


(St)o 


(BI 5 ) 
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It is Interesting to compare equation (BI5) with the empirical formula 
deduced from reference 23 « The boundary -layer thickness of reference 23? 
with u/ug 0-99 j can he rewritten 


/ 


M-s 

2pg(du/dx)< 


3.0 


(BI6) 


Equation (BI5) becomes identical with equation (B16) if the Prandtl number is 
a constant near 0.7^, which is the case treated in reference 23. 


Letting the characteristic (displacement) thickness be about one- 
half the boundary-layer thickness Sq defined by equation (BI5), one can 
define the critical Stanton number where the shock layer is dominated by the 
viscous flow and the vorticity effect so that heat conduction behind the 
shock wave and the terms in the flow equations of order (s/e) are no longer 
negligible. Thus, this critical Stanton number (St),^j^ can be given as 


(St)er 




0.62 (Xw*)* 


(BI7) 


after incorporating equation (BI5) with (6'^/Lo) I/2 (6/Lq) 1.0. 


The Blowing Parameter for the Boundary -Layer Blowoff, B^ 


As the blowing rate Increases, the effect of the foreign gases on heat 
transfer can become much stronger than that of the free -stream gas. The 
onset of this behavior is defined by "an effective blowoff parameter B*," 
described below. The relation 5 /Lo ~ L/Lq - 1 is used to replace the orig- 
inal boundary-layer thickness (d/L)^ by the increase of the shock-layer 
thickness due to the injection. Then one obtains the result from equations 
(B11) and (BI5) that 




1.82 ^ 


(BI8) 


For air as the free -stream 


gas, ( xj)* 


1/2, e — 1/16, and it follows that 


B* - 1.59 



(BI9) 


As noted from equation (BI9) , the parameter B* is independent of Stanton 
number, and thus of the product PqoR. 
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Equations for the Heat -Transfer Parameter 

The heat -transfer rate at the stagnation point for the no injection case 
can, from equation (l4), he written as 


<lco = ^ 


2 / * V 

It ' ^ w' 


J 




Pr*Re* 


(B20) 


where 


Pr* = 


(k/cp)* 


Re* = 


P V 
*^00 00 O 


M- 


and 


Poo^co ~ 


I - 1 j2(Hs - h^) 






(B21) 


(see ref. 24 for values of l/e) . It can he shown that (eqs. (l4), (B3), and 
(B21) ) 


CO JP. 


(^7/jt)(l + l/e)^^^(k/c )* 

(hg - h^) =“ const, (hg - h^) 

(h - h 

' fVl ' 


(B22) 


since (k/cp)* «:^/hs as shown in figure 2(h). Equation (B20) in dimensionless 
form is 


(St)^jR?r= JJ 


(X ')*/Pr* 


(B23) 


The present analysis on heat transfer may he expressed in term of the 
following simple correlation formula 


where 


Nu 
n/ Re 


X 


1 Wu 

\/ka Re ^ 

Ip 



X 


•Ico^^pw 

kw(hs - IV) 7 Pw^sx/p^ 


(B24) 


(B25) 


Nu 





hw(hs - h-;^) Pvj-VsLo, 


(B 26 ) 
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It follows from eq^uations (B20) and {'S)2k) that 

Nu 

X 




[(x;k/cp)^]/[(k/cp)^] 

^a(pw/Ps) 


(B27) 
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TABKE I*.- (k/cp)/(k/cp)^gf and (k/Cp)*/(k/Cp)pef 
[k/cp = (a/Pr; (k/cp)pef “ 36.^10"® Ibin/sec-ft, or = 550x10"® gm/sec-ctn] 


VsP, atm 
T, 

100.0 

10.0 

1.0 

0.1 

0.01 

1,000 

0 . lOOOOOE+Ol 

.605691E+00 

O.lOOOOOE+01 

.60569IE+OO 

O.lOOOOOE+01 

.605691E+00 

0 . lOOOOOE+01 
.605691E+00 

0 . lOOOOOE+01 
.605691E+00 

2,000 

. 145687E+01 

•923954E+00 

. I45687E+OI 
•923954E+00 

.145687E+01 

.92395^E+00 

.147018E+01 

.925063E+00 

•1555^7E+01 
• 932170E+00 

3,000 

.189730E+01 

.117341E+01 

. 20647 lE+01 
. 118814E+01 

.223923E+01 

.122759E+01 

.220755E+01 
. 1259^4E+01 

.189730E+01 

.127607E+01 

4,000 

.256478E+01 
. 144087E+01 

.249190E+01 

.148419E+01 

.218206E+01 
. 149119E+01 

. 21993 lE+01 
.147188E+01 

.270202E+01 

.149702E+01 

5,000 

.2675^6 e +01 

.168588E+01 

.2i^8227E+01 

.167247E+01 

.313709E+01 

.168791E+01 

.332438E+01 

.175475E+01 

.316304E+01 

.181001E+01 

6,000 

. 27807 OE+01 
. 18540 lE+01 

•364194E+01 
. 192837E+01 

.365864E+01 

.198598E+01 

.315177E+01 
. 20203 3E+01 

.256074E+01 
. 198097E+01 

7,000 

. 403066E+01 
.210410E+01 

.397903E+01 

.220394E+01 

.333086E+01 

.220841E+01 

.288045E+01 

.214764E+01 

.348265E+01 

.211008E+01 

8,000 

.43I926E+OI 

.236908E+OI 

.368086E+01 

.241156E+01 

.319855E+01 

.233058E+01 

. 495698E+01 
.230286E+01 

.815476E+01 

.268611E+01 

9,000 

.422441E+01 

.258319E+OI 

.35459i*'E+01 

.254053E+OI 

. 439855E+01 

.247592E+01 

.910312E+01 

.291219E+01 

.940555E+01 

.337641E+01 

10,000 

. 404048E+01 
.273753E+OI 

. 44281IE+OI 
.266996E+OI 

. 88597 lE+Ol 
.294323E+01 

.103928E+01 

.359920E+01 

.942464E+01 

.398819E+01 

11,000 

.424119E+01 

.286I78E+OI 

. 89507 8E+01 

.309010E+01 

. 104223E+02 
.355433J1+01 

. 106977E+02 

.423596E+01 

.850030E+01 

.444249E+01 

12,000 

. 491808E+OI 
.3OOO33E+OI 

. 106360E+02 
.365045E+01 

. 113205E+02 
. 416510E+01 

.102679E+02 

.476105E+01 

.764167E+01 

.474184E+01 

13,000 

.90)^886E+01 

.327562E+OI 

.118311E+02 

.423563E+01 

. 114960E+02 
.472790E+01 

.941163E+01 

.515283E+01 

. 69006 5E+01 
. 49357 lE+01 

14,000 

. II8427E+O2 
.382166E+OI 

. 127 148E+02 
.481559E+01 

.111713E+02 
. 520327E+01 

.883452E+01 
• 5^3S65E+01 

•605537E+01 

.504802E+01 

15,000 

.I31548E+O2 
. 44087 4E+01 

.131431E+02 

.535472E+01 

.105751E+02 

.557748E+01 

.836427E+01 

•564653E+01 

. 400780E+01 
.506551E+01 


Thermo -transport property of reference 49. Some changes in transport 
properties between references 3P ^9 have been noticed. Wo significant 

effects on the averaged quantity and heat -transfer parameter (table II) are 
noted. 
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TABLE II.- ^ 


atm 
T, °K^N. 

100.0 

10.0 

1.0 

0.1 

0.01 

0 

0 

0 

0 

0 

0 

1,000 

.124758 

.124758 

. 124758 

.124758 

. 124758 

2,000 

.260934 

.260934 

.260934 

.266556 

.274306 

3,000 

.424226 

.441701 

. 499412 

.587114 

.719527 

4,000 

.658638 

.793711 

.927362 

..990977 

1.04096 

5,000 

1.00186 

1.13792 

1.22421 

1.38928 

1^44939 

6,000 

1.31813 

1.49654 

1.82692 

2.62735 

3.71610 

7,000 

1.70484 

2.15229 

3^11495 

4.08941 

4.38170 

8,000 

2.29662 

3.26810 

4.30187 

4.65620 

5.24992 

9,000 

3. .18093 

4.34950 

4.87940 

5^57 992 

6.70618 

10,000 

4.17155 

5^01919 

5.66098 

6.91708 

9.29398 

11,000 

4.97112 

5.77139 

6.73418 

8.90170 

13^0915 

12,000 

5.55558 

6 . 66062 

8^5515 

11.7154 

16.7196 

13,000 

6^15677 

7 .72790 

10.0512 

15.0097 

18.8342 

14,000 

7.01739 

9.03096 

12.4241 

17.8402 

19.6352 

15,000 

7^96781 

10.5858 

15.1278 

19.6961 

19.8450 
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